Optimizing the Majorana character of SQUIDs with topologically non-trivial barriers 
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We have modeled SQUIDs with topologically non-trivial superconducting junctions and performed an op- 
timization study on the Majorana fermion detection. We find that the SQUID parameters fit,, and fie can be 
used to increase the ratio of Majorana tunneling to standard Cooper pair tunneling by more than two orders 
of magnitude. Most importantly, we show that dc SQUIDs including topologically trivial components can still 
host strong signatures of the Majorana fermion. This paves the way towards the experimental verification of the 
theoretically predicted Majorana fermion. 



I. INTRODUCTION 



Superconducting junctions with topologically non-trivial 
barriers are predicted to host Majorana bound states-^. Non- 
trivial states include the edge or surface of the recently discov- 
ered topological insulators 3-12 and semiconducting nanowires 
in the presence of Rashba spin orbit coupling and a Zee- 
man fiel d 13 ' 14 . Candidates with high potential for the de- 
tection and manipulation of the Majorana fermion-^ are su- 
perconducting quantum interference devices (SQUIDs) 16 ' 17 . 
The appearance of Majorana bound states in superconduct- 
ing junctions enables tunneling of quasiparticles with charge 
e across the junction, which doubles the Josephson periodic- 
ity, I c — Jo sin(0/2)i. The doubled periodicity is predicted 
to lead to the absence of odd integer Shapiro steps in individ- 
ual junctions, and a SQUID modulation period of 2$o instead 
of the usual $o periodicity, with <f>n = the magnetic flux 
quantum in superconductivity 1 £^£. 



Experimental efforts have been made to contact supercon- 
ductors to topologically non-trivial states^2^. Josephson ef- 
fects have been observed, and SQUIDs have been reported 22 . 
The first signatures of a Majorana fermion, characterized by a 
zero bias conductance peak, have been observed in supercon- 
ductor - semiconducting nanowire junctions 23 . Nonetheless, 
so far only $o periodic dependences have been observed. Re- 
laxation to equilibrium state a 16 ' 18 , quantum phase slips—, and 
the large bulk shunt present in contacts with topological insu- 
lators so far, may reduce the 2$o periodicity. 

A key question therefore is how to optimize the Majorana 
character. Here, we study extrinsic parameters that can be 
controlled to optimize the sin(0/2) signal from the Majorana 
fermion in dc SQUIDs composed of junctions containing both 
sin(0) and sin(0/2) components in different proportions. Our 
main observation is that the SQUID parameters and [3c are 
important parameters altering the periodicity. Furthermore, 
a superconducting interferometer will have the periodicity of 
the component with the smallest periodicity. Nonetheless, 
even in dc SQUIDs with topologically trivial components the 
Majorana character strongly influences the dc SQUID charac- 
teristics. This study is also of relevance for dc SQUIDs com- 
posed of junctions with higher order periodicities, occurring 
in SNS and SFS systems^ 4 .. 
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FIG. 1: Schematic representation of the considered dc SQUID. The 
dc SQUID is composed of a superconducting ring interrupted by 
two Josephson junctions. Charge transport through the Josephson 
junction is characterized by standard Cooper pair tunneling (sin(0)) 
and single electron tunneling by virtue of the Majorana fermion 
(sin(0/2)). The relative contribution of these two processes is de- 
termined by the factor a. We analyze two scenarios, in (a) the super- 
conductor is macroscopic and is in the trivial state, 7 P ; = 1 and in (b) 
the superconductor is in the topologically non-trivial state, y p i = 2. 
Doubled periodicity is only observed when there is a topologically 
non-trivial state in the entire ring, as in (b). However, the unusual 
current phase relation changes the dc SQUID characteristics even 
when the ring is not entirely topologically non-trivial. 



II. MODEL 



In the next session we introduce fluxoid quantization in 
superconducting rings composed of topologically trivial and 
non-trivial parts. Then we will use the fluxoid quantization 
conditions to determine the critical current in the SQUID un- 
der applied magnetic fields and derive the voltage state assum- 
ing that the junctions can be described with the resistively and 
capacitively shunted junction model. 
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A. Fluxoid quantization in topologically (non)-trivial rings 



The fluxoid quantization in a superconducting loop T leads 
to 7 P /$o periodicity, with j p i related to the charge carrier 
q = — in the loop. In macroscopic systems, j P i = 1, 

7pZ * 

but for mesoscopic systems on the order of the supercon- 
ducting coherence length £ 25 - 26 and systems including Majo- 
rana fermions, j p i can be either 1 or 2 depending on parity 
conservation 27 . Integrating the phase of a superconducting 
loop containing N Josephson junctions results in: 
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Here, T' denotes the contour of the superconducting ring 
with the Josephson junctions excluded, and A is a normaliza- 
tion constant for the current. The phase-drop over junction 



i is given by with j p j connected to the charge carrier 



q 



— in the junction. We will consider scenarios where 



7pj 
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the superconductor is either trivial 7„; = 1 or topologically 
non-trivial -f p i = 2. Also we will consider the junctions to 
be trivial 7 PJ = 1, topologically non-trivial j p j = 2 or that 
both charge carrier types are present in the junctions. When 
the junctions are topologically non-trivial, but the supercon- 
ductor is macroscopic, quantum phase slips can occur in the 
superconductor so that 7 P j can be different from j p j. Con- 
tour integration over the magnetic vector potential A results 
in the total flux Then, in the limit J s = 0, assuming 
thick superconducting leads, fluxoid quantization reduces to: 

<t- Y*, N i — H %r- = n. The flux $ is the sum of the exter- 

nal flux and the self-flux induced by current flowing through 
the ring. 

Now we will consider the case of a ring containing two 
junctions, as depicted in Fig. Q] We consider the dc SQUID 
to be symmetric, except for the current phase relationship of 
the individual junctions. Inclusion of asymmetry (e.g. in- 
ductance, critical current or capacitance asymmetry) is easily 
included. However, inclusion will only lead to asymmetri- 
cal SQUID characteristics, and will not change the period- 
icity. The total flux of the considered system is given by 
<I> = $ e + LIcXi — LI C X2- Here, $ e is the externally applied 
flux, L the inductance of a single arm and I c the critical cur- 
rent of the individual junctions 1 and 2, I c = 7 cl = I c2 . The 
factors xi,2 denote the current dependence on the phase dif- 
ference of the individual junctions, which we limit to XX® 1 = 
a sin(0) + (1 — a) sin(</>/2), with a G [0, 1] the relative am- 
plitude and xo a normalization factor to have max(x) = 1. 
The sm(4>) component is the standard Josephson relation, and 
SNS-junctions are well described by this sinusoidal relation, 
but it can include higher order components due to n Cooper 
pairs tunneling 24 , described by I s {4>) = S^Li ^" sm ( n 0n)- 
Our simplification includes the lowest frequency, which is 
enough for our conclusions. The sin(0/2) component is due 
to single electron tunneling by virtue of the Majorana fermion 
resulting in the 4tt current phase relationship periodicity. 



B. SQUID characteristics in the superconducting and voltage 

state 



The critical current for an applied external field is obtained 
by finding the solution of the fluxoid quantization equation 
with the maximal critical current. The junctions in the volt- 
age state are modeled with the resistively and capacitively 
shunted junction (RCSJ) model, assuming an ideal Joseph- 
son junction shunted by a resistor R and a capacitor C: 



I = C 



dV 
dt 
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The voltage is related to the time 



derivative of the phase by V — -§i^, the same for a topo- 
logically trivial and non-trivial ring since we have written the 
phase as - - . The RCSJ model leads to the expression 

d 2 <j> j ]_d± _ | „,,2/ 
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-w p (xi,2 — j~) = 0, with the plasma frequency 



' \IW%' ^ e SQUID parameters are defined as 
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Applying the fluxoid quantization equation, the voltage 
state can be described by the two differential equations which 
we have solved numerically 



III. RESULTS 



0. 



In this section we show the SQUID characteristics. We will 
start analyzing a dc SQUID composed of a topologically triv- 
ial ring, and two non-trivial junctions. In this regime there is 
no doubled fluxoid quantization, however the unusual current 
phase relationship causes a deviation from standard SQUID 
characteristics. After that, we will consider the SQUID in 
the entirely non-trivial regime, where doubled periodicity is 
observed due to the appearance of the Majorana fermion. Fi- 
nally, we move to the voltage state and consider both cases in 
this regime. 



A. dc SQUIDs composed of trivial and non-trivial elements 

The considered dc SQUID is composed of a topologically 
trivial ring, and two non-trivial junctions, as shown in Fig. 
QJ. Figure [2] shows the critical current dependence of the dc 
SQUID. In this figure, the non-trivial junctions develop their 
current phase relationship from pure sin(</>) to pure sin(</>/2) 
in steps 5a = 0.2. Note that all graphs are calculated for 
Pl = 0, a situation which for standard SQUIDs leads to a 
complete critical current modulation. The 2<£>o periodicity 
due to the sin(</>/2) component, tends to be completely ob- 
scured by the trivial superconducting ring. Quantum phase 
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FIG. 2: 1$ characteristics of a dc SQUID composed of a topolog- 
ically trivial superconducting ring and non-trivial junctions. The 
current phase relationship of the junctions is shown for asin(0) + 
(1 — a) sin(0/2) in steps 5a = 0.2, Pl = 0. The increase of the 
sin(<j!>/2) component causes a decrease of the oscillation amplitude 
without introducing a 2<E>o component. 



slips cause the usual $o periodicity, equivalent to what is cal- 
culated by Heck et al. 2& when one of the junctions is topo- 
logically trivial. When \i — sin(</>/2) and \2 = sin(</>) we 
obtain the result of Fig. |2]for a approximately 0.7, shifted by 
an additional j&o- Instead of 2$o periodicity, the sm(<fi/2) 
component influences the modulation depth of the SQUID. 
There is no appearance of asymmetry as is the case for asym- 
metric SQUIDs, with different critical currents of the individ- 
ual junctions 29 . The decrease in modulation depth by increas- 
ing the sin(0/2) component looks similar to increasing /?£ in 
standard SQUIDs. However, this is a parameter that can be 
controlled externally, and a large f3 L results in more triangular 
oscillations. If one junction is topologically trivial, the same 
effect occurs, combined with a phase shift due to asymme- 
try between the junctions. Therefore, even in rings including 
topologically trivial components, a sin(</>/2) current phase re- 
lationship can be detected, although the effect is more subtle 
than a 2$o periodicity. 
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FIG. 3: 1$ characteristics of a dc SQUID composed of a topologi- 
cally non-trivial superconducting ring and non-trivial junctions, (a) 
dc SQUID oscillations for two symmetric junctions with equal am- 
plitude sin(0) and sin(0/2) components, a — 0.5. Increasing Pl 
(in steps SPl = 0.5, and shifted for clarity) promotes the 2$o pe- 
riod, since the effective screening is smaller for Majorana tunneling 
than Cooper pair tunneling, (b) (blue) FFT amplitude of the two 
components as function of Pl- (red) Evolution of the FFT amplitude 
ratio as a function of Pl- The dashed lines represent the ratio when 
the junctions have only 5% sin(</)/2) component contribution. The 
2<I>o component can be more than 2 orders of magnitude larger than 
the <J>o component. 



B. Topologically non-trivial SQUIDs 

In the case when the ring is completely topologically non- 
trivial, corresponding to Fig. QJ), 2$ periodicity is to be ob- 
served. In the limiting case [3 L = and / = I c sin(0/2), 
the critical current dependence on field can be written as 
/ = 2I C \ cos(7T2§-)|, resulting in the 2$ periodicity. If both 
sin(0) and sin(0/2) components are present in the junctions, 
both periodicities are observed, as shown in Fig. [3^ for equal 
ratios in the junctions. Interestingly, increasing f3 L results in 
a larger 2$o component and a reduced $o component. In Fig. 

the ratio dependence on f3 L is shown, where the ratio is de- 
fined using the frequency amplitude after Fourier transforma- 



tion. The screening parameter (3 L is composed of the critical 
current of the junctions, and the inductance of the ring de- 
termined by geometrical factors, but is also dependent on the 
charge carrier in the ring. Since Majorana tunneling is with 
charge e instead of 2e as is the case for Cooper pair transport, 
the effective screening is reduced by a factor 2. By optimizing 
f3 L using the tunable inductance of the ring, it is therefore pos- 
sible to dramatically increase the 2<E>o component relative to 
the standard $o-P er i°dic component, ideal for the observation 
of the Majorana fermion. 
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FIG. 4: Topologically non-trivial dc SQUID in the voltage mode, (a) dc SQUID IV characteristics for junctions with sin(<jf>) and sin(</>/2) 
current phase relationships; fie = l,fit=l. The sin(</>/2) (solid line) component doubles the resonance voltage V res (dashed line) with respect 
to the standard SQUID. The IV characteristic of the standard SQUID is shifted for clarity, (b) Damping and resonance in the SQUID. The 
sin(<£/2) component oscillates with half the frequency, so that the sin(^) and sin(</>/2) components have their minima and maxima at different 
voltages (solid red are maxima for the sin(0/2) component and dashed blue maxima for the sin(<£) component, (c) V<f> characteristics of a dc 
SQUID with equal sin(0) and sin(0/2) components in the junctions, and /3c = 1, Pl = 1. The SQUID shows both $0 and 2$o periodicity, 
depending on the voltage close to V r e S ,4> or V reStt j,/2, which is controlled by the bias current. 



C. The voltage state 

When the SQUID is operated in the voltage mode, both $0 
and 2$o periodicity can be observed, even if the ring includes 
trivial components. The relative amplitude depends on the 
voltage, controlled by the bias current. Figure |4fc shows the 
IV characteristics for a dc SQUID with either pure sin(^) or 
sin (0/2) components. The current modulation for high volt- 
ages is inverted with respect to the modulation for small volt- 
ages. This is the result of the nonlinear interaction of the ac 
Josephson current with the resonant circuit formed by the loop 
inductance L and the junction capacitance C. This resonance 

voltage is given by V res = 7pj \//3^r^ c ^' ^ s a resu ^' tne 
sm((f>) and sin(0/2) components cause oscillations with dif- 
ferent frequency as a function of voltage, see Fig. |4j3. Conse- 
quently, the relative amplitude of the 2 and 4tt periodicity of 
the V((f>) characteristics, shown in Fig. @J for pc — 1> an d 
Pl = 1 depend on the resonance voltages V re s,4> and V res ^/2, 
which is controlled by the bias current. The resonance voltage 
is independent of the loop parity, but for 7 p ; = 1 the oscilla- 
tion amplitude is reduced by quantum phase slips. This is sim- 
ilar to the cause of the incomplete critical current modulation 
shown in Fig. [2] 

The damping of the current modulation with increasing bias 

voltage is characterized by Q — J independent on the 

current phase periodicity, and SQUIDs with sin(0/2) com- 
ponent junctions have therefore a smaller current modula- 
tion where the modulation is inverted, compared to standard 
SQUIDs. Choosing a large /3l will decrease the damping term 
and increase the amplitude of the voltage resonances. A small 
f3c will reduce the damping term but increase the resonance 
voltage. 



IV. APPLICATIONS TO TOPOLOGICALLY 
NON-TRIVIAL SYSTEMS 

We now discuss the implications of our proposals. The pos- 
sible realization of a new emergent particle in condensed mat- 
ter physics together with the potential for quantum compu- 
tation, boosted the search for superconducting systems host- 
ing the Majorana fermion. Superconductor - semiconductor 
structures in the presence of strong spin orbit coupling and 
Zeeman fields are currently quite successful and signatures 
of Majorana fermions have been observed characterized by 
zero bias conductance peaks 23 . In superconductor - topo- 
logical insulator junctions, Josephson supercurrents have also 
been observe d 19-21 ' 30 . In the case of topological insulator sys- 
tems, bulk shunting likely introduces sin(</>) terms in the cur- 
rent phase relation. Quantum phase slips 17 and quasiparticle 
poisoning 16 will be relevant for all proposals, relaxing the sys- 
tem to sin(0) periodicity. Therefore, increasing the sin(^>/2) 
component is important for all these proposals. 

The parameter /3l is easily tunable. This SQUID parameter 
depends on the inductance, determined by the SQUID geom- 
etry, and independent of the individual junctions. The param- 
eter /3c is a junction parameter and therefore more difficult to 
tune. However, varying the length and width of the junctions 
and controlling the interface transparency tune /3c- 



V. CONCLUSIONS 

In conclusion, we have studied dc SQUIDs containing two 
topologically non-trivial barriers. The 2<£>o periodicity stem- 
ming from the Majorana fermion can only be detected in non- 
trivial dc SQUIDs. However, even loops containing topo- 
logically trivial elements are influenced by the presence of 
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junctions with sm(4>/2) components. This is observed both 
in the critical current modulation by flux and the resonance 
voltage. The SQUID parameters can be used to increase the 
relative component; increasing /3j, is found to largely increase 
the component with the largest periodicity. The V{<f) rela- 
tion is altered when both components are present, and both 
components can be maximized at different bias currents, de- 
termined by the resonance voltage. In recently fabricated S- 
TI-S junctions, fie is usually low due to bulk shunting. This 
increases the resonance voltage to a regime where the damp- 
ing is higher, which complicated observing a clear differ- 
ence between <£>o an d 2$o periodicities. Decreasing this bulk 
shunt would therefore simplify the 2$o periodicity observa- 
tion. Nonetheless, the strong effect of fit,; large (3l increases 



the 2$o component over 100 times, allows the detection of 
the Majorana fermion, under the right intrinsic circumstances 
(e.g. no relaxation to equilibrium), in dc SQUIDs composed 
of present day S-TI-S junctions. This result is also of rele- 
vance to devices where the Majorana character is induced via 
other means, such as in semiconducting nanowires with strong 
Rashba spin orbit coupling. Tuning /3l in combination with 
ac measurements to prevent relaxation paves the way to the 
observation of exotic properties of the Majorana fermion. 
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